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Pressure and higher-order spectra 
for homogeneous isotropic turbulence 

By D. I. Pullin 1 AND R. S. Rogallo 2 


The spectra of the pressure, and other higher-order quantities including the dis- 
sipation, the enstrophy, and the square of the longitudinal velocity derivative are 
computed using data obtained from direct numerical simulation of homogeneous 
isotropic turbulence at Taylor- Reynolds numbers R\ in the range 38 - 170. For the 
pressure spectra we find reasonable collapse in the dissipation range (of the velocity 
spectrum) when scaled in Kolmogorov variables and some evidence, which is not 
conclusive, for the existence of a k~ 7 / 3 inertial range, where k = |k| is the modu- 
lus of the wavenumber. The power spectra of the dissipation, the enstrophy, and 
the square of the longitudinal velocity derivative separate in the dissipation range, 
but appear to converge together in the short inertial range of the simulations. A 
least-squares curve-fit in the dissipation range for one value of R\ = 96 gives a form 
for the spectrum of the dissipation as k° exp(— Cfcr?), for kt) > 0.2, where r) is the 
Kolmogorov length and C « 2.5. 


1. Introduction 

The collapse of data for the various forms of the velocity power spectrum when 
scaled in Kolmogorov variables of the form 

E(k) = (^ 5 { £ )) 1/4 /i (*«?), (1) 

is well established as either an exact, or at least a very good approximation at 
large Reynolds numbers, by both experiment and numerical simulation; see for 
example Saddoughi & Veeravalli (1994) for an update of a graphical compilation 
originally due to Chapman (1979). In (1) E is either the shell-summed spectrum or 
a form of the related one- dimensional spectra, r) = (i/ 3 /^)) 1 / 4 is the Kolmogorov 
length, and (e) is the volume- averaged dissipation. The aim of the present work 
is to use numerical data bases obtained from direct numerical simulation (DNS) of 
homogeneous isotropic turbulence to study the power spectra of several quantities 
that are quadratic in the velocity or its spatial derivative. It is hoped that the 
results may be useful for testing the predictive capability of theories or models of 
turbulent fine scales as well as adding to our basic understanding of turbulence in 
the inertial and dissipation ranges. 
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We shall consider two-point correlations and associated power spectra of quanti- 
ties that are scalar functions of position in homogeneous isotropic turbulence. Let 
q(Xjt) be an arbitrary function of position x = (x\,X 2 )X$) and time t, and denote 
its Fourier transform by 

m = J Jq(x)e- k "dx. (2) 

In (2) and subsequently we have suppressed the explicit dependence on t. Define 
the two-point one-time correlation of <?(x) at points x, x + r by 

Q( r) = («(x)g(x + r)), (3) 


where (...) is a volume average over x. The Fourier transform of Q(r) is 

n(k) = 5 (4) 


and it is easy to show that 

n(k) = 9(k)9*(k), (5) 

where * denotes the complex conjugate. On integrating (5) over the surface of a 
sphere S(k) of radius k in k-space, we obtain 


E q (k) = 


I Is(k) 


n{k)ds(k). 


( 6 ) 


We shall refer to g(k) q*( k) as the power spectrum of q(x) and to E q (k) defined by 
(6) as its shell-summed power spectrum. When the turbulence is isotropic and <?(x) 
is an invariant scalar in the sense that its value at any point is axis independent, 
then Q( r) = Q(r), r = |r|, and II(k) = II(fc), k = |k|. The right-hand side of (6) is 
then equal to 47rfc 2 II(fc). 

We consider power spectra of four quantities for which <?(x) is identified in turn 
with the following: the pressure-density ratio p(x), the dissipation e = 2u(du t / dxj + 
duj/dxi ) 2 , the enstrophy density ft = a> 2 , and a quantity proportional to the square 


of the longitudinal velocity derivative T = 15 (du) 2 , where du = du\jdx\ and 
u = (ui,ii 2 ,U 3 ) is the velocity vector. We denote the shell-summed power spectra 
of these quantities, in Kolmogorov scaling variables, respectively by 

E p (k) = {ef*u 7 '*Mkri), (7) 

E e (k ) = (e) 7/i u 3/i f 3 (kTj), (8) 

E n (k) = ( £ ) 7/ V- 5 / 4 /4(*»?), (9) 

£ T (fc) = (e) 7/ \-^ 4 Mkv)- (10) 
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For homogeneous isotropic turbulence, (e), (fl), and (T) are nonzero and satisfy 
the relations 

(e) = i/(0) =i/( T). (11) 

In (8-10) these mean values axe subtracted, and the integrals of the right-hand sides 
over k = (0, oo) are equal to the fluctuations, respectively, = ((e — {e}) 2 }, 

(O' 2 ) = ((0 - (0)) 2 ), and (r 2 ) = ((T - (T)) 2 ). 

We remark that (du\ / dxi) 2 is not a proper scalar, but is one component of the 
fourth-order tensor 

\dxj) \dx m J 

whose two-point correlation is an eighth-order tensor. As a consequence the power 
spectrum of T is not uniform over spherical surfaces in k-space. We nevertheless 
treat T as if it were a scalar for the reason that it has been used as a surrogate for 
€ in experiment, being the gradient-variance that can most easily be measured with 
a hot-wire probe. It is therefore of interest to compare its fluctuations with those 
of e and fi. Note that (3-5) remain valid and (6) may still be applied. 

2. Data bases 

For the most part we have used data from the forced DNS runs of Jimenez et al 
(1993), henceforth referred to as JWSR. See their §2 for a detailed description of 
the numerical experiments and their Table 1 for a summary of the flow parameters. 
Our data base differs from JWSR Table 1 as follows: first, JWSR take both volume 
and time averages, the latter being over several large-eddy turnover times, whereas 
our results are based on Fourier coefficients of the velocity from a single time frame, 
being the last for the run. Comparisons made with one-time results at earlier 
times separated by a substantial fraction of a large-eddy turnover time indicated 
that statistical equilibrium had been obtained. Because we are using a single time 
frame our values R x = 38(64 3 ), 65(128 3 ), 96(256 3 ), and 170(256 3 ) (the bracketed 
number gives the numerical grid resolution N 3 ) differ slightly from those of JWSR. 
Secondly and more importantly, our largest R x run was at a resolution of 256 3 
compared to the 512 3 run of JWSR which was unavailable to us. Since the 256 3 
data at R x = 170 has k max r] % 1.0 compared with k max r) « 2.0 for all the JWSR 
runs, it may be somewhat under-resolved. All shell-summed spectra were calculated 
directly from coefficients in k-space with de-aliased quadratic products evaluated in 
physical space. Those for the pressure were obtained from a solution of the relevant 
Poisson equation. 

For the pressure spectra, in addition to the DNS runs, we also consider results 
from two kinds of LES (large-eddy simulation) runs. The first of these utilizes a 
method in which, at every time step, the magnitudes of the Fourier components of 
the velocity are adjusted, with no change to either the relative magnitudes in the 
(xi,X 2 ,x 3 ) directions or the complex phases, such that the velocity spectrum was 
locked to a k ~ 5 ^ 3 form. We will refer to this as Locked-E LES (She & Jackson, 
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170, Rx = 96, Rx = 65, Rx = 38. 


1993). In addition some pressure spectra calculated from a Smagorinsky subgrid 
approach using the so called dynamic localization model ( DLM , Ghosal et al. 1994) 
Eire also presented. 

3. Spectrum of the pressure 

Fig. 1 shows the velocity (energy) spectrum plotted in the compensated 
inertial-range form (e)~ 2/3 k 5 / 3 E{k) = {ki]) 5 / 3 fi(krj). The data, including the 
R x = 170 case, collapses reasonably across the whole range of kr). A Kolmogorov 
constant of £q as 2.2 appears to be indicated, but as pointed out by Jimenez (pri- 
vate communication), the plateau at kr\ < 0.2 may in fact be part of a bump in the 
spectrum near the beginning of the dissipation range (see Saddoughi & Veeravalli, 
1994, for related experimental evidence), which is not properly resolved owing to the 
small inertial range of the DNS. The above quoted K o may then be an overestimate. 

It is well known that Kolmogorov-type dimensional arguments suggest a form for 
E p in the inertial range 

E p (k) = K p (e) 4/3 AT 7 ' 3 , (12) 

where IC P is a dimensionless number which may or may not be a universail con- 
stant. In Fig. 2 E p (k) is plotted in the compensated form suggested by (12), 
(e)~ 4/3 i <c 7/3 E p (k ) = (krj) 7/3 / 2 (l-r?). The collapse in the dissipation range (by which 
we mean kr] greater than about 0.125) is only fair if one discards the Rx = 170 
case as under-resolved, and poor if it is included. If this case is disceirded Fig. 2 is 
consistent with a possible plateau corresponding to IC P « 8.5, which is somewhat 
higher than the value AC P « 7 suggested by Pumir (1994). Again, owing to the small 
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FIGURE 2. Compensated pressure spectra in Kolmogorov units. R\ = 170, 

R x = 96, R\ = 65, R\ = 38. 

inertial range, definitive conclusions cannot be drawn and we cannot rule out the 
possibility of a bump analogous to that found found for E(k) with or without a 
plateau at lower hr). 

An argument due to Obukov (1949) and Batchelor (1951) based on a joint-normal 
hypothesis for the two-point probability distribution of the velocity field suggests 
that Kq and fC p (if they exist) are related as 

= 4 ^| ICl = 1-3245105. Xjj. (13) 

Taking Kq = 1.5 gives, from (13), K p = 2.98 which is much lower than the peak 
shown in Fig. 2 but which agrees approximately with the range of values 2.4 — 3.4 
estimated by Pullin (1994) from the Lundgren (1982) stretched spiral vortex model. 

The ratio E p (k)/(k E(k) 2 ) is shown in Figs. 3-4. If both E(k) and E p (k) follow 
their respective K-41 power laws in the inertial range, this ratio should plateau at 
K p /Kq. As a test of our means of calculating E p , the phases of the Fourier co- 
efficients of the velocity field for the Locked- E LES run were randomized without 
change to their magnitudes. This should produce a Gaussian joint-normal veloc- 
ity probability distribution while maintaining a k ~ 5 / 3 inertial range. The plateau 
shown in Fig. 3 for this case is in good agreement with (13) as indeed it should be 
(see the article by Gotoh &; Rogallo in these proceedings for a discussion of various 
spectra for Gaussian velocity pdf’s). The roll-off at both large and small k is at- 
tributed to the finite range of k in the calculation of the pressure via the (elliptic) 
Poisson equation, which evidently “feels” the bounded dimensions of the box. The 
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FIGURE 5. Compensated one-dimensional pressure spectra in Kolmogorov units. 

R x = 170, R x = 96, R x = 65, Rx = 38. x DLM LES 

(64 3 ), Rx = 269. + DLM LES (32 3 ), Rx = 176. • experiment of George et al 
(1984). 

true Locked-E LES also shows a plateau in Figs. 3-4, but at a value close to 3.8. 
In Fig. 4 the DNS and Locked-E LES show agreement in a small range of k at the 
larger Rx . 

An attempt to bring experiment, DNS, and LES results together is made in Fig. 5 
which plots the one-dimensional form of the pressure spectrum 

£p 1} (*i) = l (14) 

L Jki K 

The measurements of George et al (1984) shown in Fig. 5 were made in the 
mixing layer of a turbulent jet at 1.5 — 3 diameters downstream of the jet exit. 
Reynolds numbers based on the exit velocity and jet diameter were 4,0 x 10 5 and 
6.2 x 10 5 . Taylor Reynolds numbers were in the range Rx & 350 — 600. George et 
al nondimensionalized their data with large-scale parameters which may be more 
appropriate at the fairly low wavenumbers of the experiment. We have rescaled 
their results using estimates of the dissipation and of the Kolmogorov microscale 
given in §17 of their paper. They state that velocity contamination caused by 
flow-probe interactions may have had some effect on the accuracy of the data. 
Fig. 5 well illustrates the difficulties of studying pressure correlations. Neither 
DNS, which is currently limited to low Rx , nor experiment, which is limited by 
accuracy problems at large wavenumbers, can span an appreciable spectral band. 
The combined experiment, DLM LES, and DNS results may indicate a plateau in 
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FIGURE 6. Power spectra of the dissipation in Kolmogorov units. Rx = 170, 

R x = 96, Rx = 65, R\ — 38. 

the fc^-compensated spectra, but the scatter is very large. If (12) is true then 
from (14) the one-dimensional inertial-range coefficient is ACj, 1 ^ = K p . 


4. Spectra of e, ft, and T 

Calculated dissipation spectra E t and enstrophy spectra Eq are shown in Figs. 6 
and 7 respectively. Fig. 8 shows a comparison of these together with the spec- 
trum of T = lb{du\ j dx\ ) 2 at R\ = 96. This indicates that the fluctuations of 
ft are rather larger than those of e/v with fluctuations of 15(dui/dxi) 2 larger 
still. If Kolmogorov scaling is exact then the quantities (e 2 )/ (e) 2 , (ft 2 )/(ft) 2 , 
and ((du) 4 ) / ((du)) 2 should each be constant independent of Reynolds number. 
There is evidence from DNS, e.g. Kerr (1985), JWSR, that this is not the case for 
((du) 4 ) / ((du)) 2 , which exhibits a weak dependence on the Taylor Reynolds number 
Rx = uX/v, where u is the root-mean square of one component of the velocity and A 
is the Taylor microscale. This is confirmed by Gotoh & Rogallo in these proceedings 
who also tabulate (e 2 )/(e) 2 , (ft 2 )/(ft) 2 for the present data set. 

It is notable that Fig. 8 shows no sign of power- law behavior at low k, and indeed 
an asymptote to a constant appears to be indicated for all three quantities. This 
is in qualitative agreement with measurements, made using a multi-wire probe, of 
Tsinober et al (1992) (see their Fig. 9(a)). In Fig. 9 we have plotted the spectra 
of Fig. 8 in the form d[log(f n (kr]))]/d[log(kT))}, n = 3 , 4,5 where /„ are defined by 
Equations (8-10). If curves are of the form (krj)P exp (—Ckrj), then this plot should 
be a straight line of slope -C that intercepts the vertical axis at /3. The fit of a 
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FIGURE 7. Power spectra of the enstrophy in Kolmogorov units. R\ = 170, 

R\ = 96, R\ = 65, R\ = 38. 



FIGURE 8. Comparison of quadratic velocity-derivative spectra at R\ = 96. 
e, ft, T = \b(dux / dxx) 2 . 
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k 7] 

Figure 9. Spectra of Fig. 8 plotted in the form d[log(/ n (fc?/))]/d[log(A:^)], n = 

3,4,5 where /„ are defined by Equations (8-10). e, ........ T = 

15(dui/dxi ) 2 . 

straight line to the data sets in Fig. 9 is very approximate. Nevertheless a least- 
squares best fit was made in each case in the window 0.2 < k < 1.8. The method 
was tested using the velocity spectrum (not shown) for which we found ft = — 1.6 
and C = 4.9 in fair agreement with Kerr (1990) and others. We find that for E € \ 
f3 « 0., C w 2.5, Eq ; ft ^ 0.5, C « 3.5. There was some sensitivity to the chosen 
range, so the errors are substantial, of order ±0.1 for /? and ±0.2 for C . Clearly a 
larger resolved dissipation range is needed to improve accuracy. 

5. Concluding remarks 

We have examined the power spectra of several quantities using numerical data 
bases from DNS and LES box turbulence, focusing mainly on flow variables that are 
quadratic in the velocity components and which have proven difficult to measure 
experimentally. When taken in conjunction with the one- dimensional spectra of 
Figure 5, the DNS shell-summed spectra provide tentative evidence for the presence 
of a Jt -7 / 3 range with A Z p ss 8.5 and a corresponding ICp 1 * ss 1.8, but with large error. 
These are higher than those predicted by either the joint-Gaussian model (AC P ~ 2.98 
corresponding to K . o = 1.5) or the stretched- spiral vortex model (AC P = 2.4 — 3.4 
depending on R\). Since it is unlikely that the spectral range of DNS will increase 
sufficiently in the near future to resolve this issue, there is clearly a need for a 
definitive experiment of the Saddoughi & Veeravalli (1994) type. The spectra of (e) 
and (Q.) show significant differences in the dissipation range but appear to behave 
similarly in the short inertial range, where a k° form is indicated, in agreement with 
the prediction of the stretched spiral-vortex model (Pullin et al. 1994). 
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